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Abstract — We consider the problem of rate and power allo- 
cation in a fading multiple-access channel. Our objective is to 
obtain rate and power allocation policies that maximize a utility 
function defined over average transmission rates. In contrast 
with the literature, which focuses on the linear case, we present 
results for general concave utility functions. We consider two 
cases. In the first case, we assume that power control is possible 
and channel statistics are known. In this case, we show that the 
optimal policies can be obtained greedily by maximizing a linear 
utility function at each channel state. In the second case, we 
assume that power control is not possible and channel statistics 
are not available. In this case, we define a greedy rate allocation 
policy and provide upper bounds on the performance difference 
between the optimal and the greedy policy. Our bounds highlight 
the dependence of the performance difference on the channel 
variations and the structure of the utility function. 



I. Introduction 

Dynamic allocation of communication resources such as 
bandwidth or transmission power is a central issue in multiple 
access channels in view of the time varying nature of the 
channel and interference effect. Most of the existing literature 
focuses on specific communication schemes such as TDMA 
(time-division multiple access) [1] and CDMA (code-division 
multiple access) [2], [3] systems. An exception is the work by 
Tse et al. [4], who consider the notion of throughput capacity 
for the fading channel with Channel State Information (CSI). 
That is the notion of Shannon capacity applied to the fading 
channel, where the codeword length can be arbitrarily long 
to average over the fading of the channel. The points on the 
boundary of the capacity region are attained by dynamically 
allocating the resources with the goal of maximizing a linear 
utility function. 

In this paper, we consider the problem of rate and power 
allocation in a multiple access channel with perfect CSI. 
Contrary to the linear case in [4], we consider maximiz- 
ing a general utility function of transmission rates over the 
throughput capacity region. Such a general concave utility 
function allows us to capture different performance metrics 
such as fairness or delay (c.f. Shenker [5], Srikant [6]). Our 
contributions can be summarized as follows. 
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We first consider the case where channel statistics are known 
and power can be controlled at the transmitters. Owing to 
strict convexity of the capacity region, we show that the 
resource allocation problem for a general concave utility is 
equivalent to another problem with a linear utility. Hence, the 
optimal resource allocation policies are obtained by applying 
the results in [4] for the linear utility. Given a general utility 
function, the conditional gradient method is used to obtain 
the corresponding linear utility. Second, we consider the case 
where the transmitters do not have the power control feature 
and channel statistics are not known. In this case, a greedy 
policy which maximizes the utility function for any given 
channel state, is suboptimal. However, we can bound the 
performance difference between the optimal and the greedy 
policies. We show that this bound is tight in the sense that 
it goes to zero either as the utility function tends to a linear 
function of the rates or as the channel variations vanish. 

Other than the papers cited above, our work is also related 
to the work of Vishwanath et al. [7] which builds on [4] and 
takes a similar approach to the resource allocation problem for 
linear utility functions. Other works address different criteria 
for resource allocation including minimizing the weighted sum 
of transmission powers [8], and considering Quality of Service 
(QoS) constraints [9]. In contrast to this literature, we consider 
the utility maximization framework for general concave utility 
functions. 

The remainder of this paper is organized as follows: In 
Section II, we introduce the model and describe the capacity 
region of a fading multiple-access channel. In Section III, we 
address the resource allocation problem with power control 
and known channel statistics. In Section IV, we consider the 
same problem without power control and channel statistics. 
Finally, we give our concluding remarks in Section V. 

Regarding the notation, we denote by xi the i-th component 
of a vector x. A vector x is positive when xi > for all 
components i of x. We denote the nonnegative orthant by 
R™, i.e., R'l = {x € R" | x > 0}. We write x' to denote the 
transpose of a vector x. 

II. System Model 

We consider M transmitters sharing the same media to 
communicate to a single receiver. We model the channel as 
a Gaussian multiple access channel with flat fading effects 

M 

Y(n) = y/Hi{n)Xi(n) + Z(n), (1) 
i=i 

where Xi(ri) and flj(n) are the transmitted waveform and the 
fading process of the i-th transmitter, respectively, and Z(n) 
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is white Gaussian noise with variance Nq. We assume that the 
fading processes of all transmitters are jointly stationary and 
ergodic, and the stationary distribution of the fading process 
has continuous density. We also assume that all the transmitters 
and the receiver have instant access to channel state infor- 
mation. In practice, the receiver measures the channels and 
feeds back the channel information to the transmitters. The 
implicit assumption in this model is that the channel variations 
are much slower than the data rate, so that the channel can 
be measured accurately at the receiver and the amount of 
feedback bits is negligible compared to that of transmitting 
information. 

First, consider the non-fading case where the channel gains 
are fixed. The capacity region of the Gaussian multiple-access 
channel with no power control is described as follows [10] 

C g (P,h) = iReRM :J2Ri<c(j2h i P i ,N 

for all S CM = {1,..., M}\, (2) 



where P, and Ri are the i-th transmitter's power and rate, 
respectively. C(P, N) denotes Shannon's formula for the ca- 
pacity of AWGN channel given by 



C(P, N) = ilog(l + ^) nats. 



(3) 



For a multiple-access channel with fading, but fixed trans- 
mission powers Pi, the throughput capacity region is given by 
averaging the instantaneous capacity regions with respect to 
the fading process [11], 



C a {P) = IreRM :^2Ri<E H cf^HtPi,^ 



for all S C M 



Cf(V) 



Re 



E 



H 



, for all S C M 



Tse et al. [4] have shown that the throughput capacity of a 
multiple access fading channel is given by 



C(P) = |J C f (V), 



(5) 



where Q = {V : E H [P Z (H)} < P 4 ,for all i} is the set of all 
power control policies satisfying the average power constraint. 
Let us define the notion of boundary or dominant face for any 
of the capacity regions defined above. 



Definition 1: The dominant face or boundary of a capacity 
region, denoted by J-(-), is defined as the set of all M-tuples 
in the capacity region such that no component can be increased 
without decreasing others while remaining in the capacity 
region. 

III. Rate Allocation with Power Control 

In this section, we assume that the channel statistics are 
known a priori. The goal of optimal resource allocation is to 
find feasible rate and power allocation policies denoted by 1Z* 
and V*, respectively, such that K*(H) € C g (P*(H),H), 
and V* G Q. Moreover, 

E H [n*(H)}= R* = argmax u(R) 

subject to Re C(P) (6) 

where u(-) is a given utility function and is assumed to 
be a continuously differentiable concave function of R, and 
monotonically increasing in each component Ri for all i. 

For the case of a linear utility function, i.e., u(R) — /jl'R 
for some /i e R+ f , Tse et al. [4] have shown that the optimal 
rate and power allocation policies are given by the optimal 
solution to a linear program, i.e., 



(4) 



where H is a random vector with the stationary distribution 
of the fading process. 

A power control policy V is a map from any given fading 
state h to P(h) = ("Pi (ft,), . . . ,VM(h)), the powers allocated 
to the transmitters. Similarly, we can define the rate allocation 
policy, 1Z, as a map from the fading state h to the transmis- 
sion rates, IZ(h). For any given power-control policy V, the 
capacity region follows from <j4j as 



(K*(h),P*(h)) = argmax (//r - \'p) 

r ,p 

subject to r e C g (h 7 p), (7) 
where h is the channel state realization, and A e TC u 



is a 



Lagrange multiplier satisfying the average power constraint, 
i.e., A is the unique solution of the following equations 



h j 2\ i (N + z) 



IP* 



2X k h(N + z) 



2Xi(N + z) + (fi k - m)h 



f l {h)dh&z = P l 



(8) 



where F k and f k are cumulative distribution function (CDF) 
and probability density function (PDF) of the stationary 
distribution of the channel state process for transmitter k, 
respectively. 

Exploiting the polymatroid structure of the capacity region, 
problem (0 can be solved by a simple greedy algorithm 
(see Lemma 3.2 of [4]). It is also shown in [4] that for 
positive /j, the optimal solution, R*, to the problem in (O 
is uniquely obtained. Given the distribution of channel state 
process, denoted by F k and f k , we have 



1 



n 



o 2(iVo + z) y gMga±g] 
2X k h(N + z) 



Fi 



2Xi(N + z) + (jj, k - m)h 



fi(h)dhdz, 



(9) 



The uniqueness of R* follows from the fact that the stationary 
distribution of the fading process has continuous density [4]. 
It is worth mentioning that (O parametrically describes the 
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boundary of the capacity region, and hence, there is a one-to- 
one correspondence between the boundary of C(P) and the 
positive vectors fi with unit norm. 

Now consider a general concave utility function. We use 
an iterative method to compute the optimal solution, R*, of 
problem ©. Note that by mono tonicity of the utility function, 
R* always lies on the boundary of the capacity region, C(P). 
Once R* is known, then in view of one-to-one correspondence 
between the boundary of C(P) and the positive vectors fi, 
there exist a positive vector fi* such that 

H*= argmax (n*)'R subject to ReC(P). (10) 

Therefore the optimal rate and power allocation policies can 
be obtained by using the greedy policies of Tse et al. [4] for 
the linear utility function, u(R) = (fi*)'R. 

We use the conditional gradient method [12] in order to 
iteratively compute the optimal solution of problem ©. The 
k-th iteration of the method is given by 



R k+l =R k + Q k( R k 

where a k is the stepsize and R k is obtained as 

R k e argmax (vu(R k )'(R - R h 
Rec(p) V 



(11) 



(12) 



Since the utility function is monotonically increasing, the 
gradient vector is always positive and, hence, the unique 
optimal solution to the above sub-problem is obtained by 
(O, in which fi is replaced by \7u(R k ). By concavity of 
the utility function and convexity of the capacity region, the 
iteration (fTTT > will converge to the optimal solution of © for 
appropriate stepsize selection rules such as Armijo rule or 
limited maximization rule (c.f. [12] pp. 220-222). 

Note that our goal is to determine rate and power allocation 
policies. Finding R* allows us to determine such policies 
by the greedy policy in (0 for fj,* — Wu(R*). It is worth 
mentioning that all the computations for obtaining R* are 
performed once in the setup of the communication session. 
So the convergence rate of the conditional gradient method is 
generally not of critical importance. 

IV. Rate Allocation without Power Control 

In this section we assume that the channel statistics are 
not known and that the transmission powers are fixed to P. 
In practice, this scenario occurs when the transmission power 
may be limited owing to environmental limitations such as 
human presence, or limitations of the hardware. 

The capacity region of the multiple access channel for this 
scenario is a polyhedron and is given by Similarly to the 
previous case, the optimal rate allocation policy, 7Z* (•), is such 
that 11* {H) e C g (P,H), and 

E H [Tl*{H)}= R* e argmax u(R) 

subject to Re C a (P). (13) 

It is worth mentioning that the approach used to find the 
optimal resource allocation policies for the previous case does 
not apply to this scenario, because C g (P, h) is a polyhedron 



and hence, the uniqueness property of R* for any positive 
vector does not hold anymore. 

Here we present a greedy rate allocation policy and com- 
pare its performance with the unknown optimal policy. The 
performance of a particular rate allocation policy is defined 
as the utility at the average rate achieved by that policy. The 
greedy policy, denoted by TZ(-), optimizes the utility function 
for any channel realization, i.e., 

K(h) = aigmax ReCs{ph) u(R). (14) 

Consider the following relations 



i H [u(K*(H))] 



< 
< 
< 



E H [u(K(H))] 
u(E H [H(H)]) 
u(E H [lZ*(H)]), 



(15) 



where the second inequality follows from the Jensen's inequal- 
ity by concavity of the utility function. 

In the case of a linear utility function we have 
u(E H [ll*{H)]) = E H [u(Tl*(H))], so equality holds 
throughout in ( TT~5T > and TZ(-) is indeed the optimal rate alloca- 
tion policy. For nonlinear utility functions, the greedy policy 
can be strictly suboptimal. 

However, the greedy policy is not arbitrarily worse than 
the optimal one. In view of (T5[ , we can bound the per- 
formance difference, u(R*) — u(Eh\TZ(H)\\ by bounding 

u(E H [TZ*(H)])-u(E H [n(H)])\ or \u(E H [R,*(H)]) - 

Eh [u (lZ*(H))] from above. We show that the first bound 
goes to zero as the channel variations become small and the 
second bound vanishes as the utility function tends to have a 
more linear structure. 

Before stating the main theorems, let us introduce some 
useful definitions and lemmas. 

Definition 2: Let Q be a polyhedron described by a set of 
linear constraints, i.e., 



{x e W 1 : Ax < b} . 



(16) 



Define the expansion of Q by 8, denoted by £g(Q), as the 
polyhedron obtained by relaxing all the constraints in (TTST l. i.e., 
£ S (Q) = {x e R" : Ax < b + SI} , where 1 is the vector of 
all ones. 

Definition 3: Let X and Y be two polyhedra described by 
a set of linear constraints. Let £d{X) be an expansion of X by 
relaxing its constraints by d. The distance dn {X, Y) between 
X and Y is defined as the minimum scalar d such that X C 
£ d (Y) and Y C £ d (X). 

Lemma Q] extends Chebychev's inequality for capacity re- 
gions. It states that the time varying capacity region does not 
deviate much from its mean with high probability. 

Lemma 1: Let H be a random vector with the stationary 
distribution of the fading process with mean H and covariance 
matrix K. Then 



Pr 



[d H (Cg(P,H),C a (P))> 5} 



< 



s 2 ■ 



(17) 
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where a 2 H is defined as 



= 1 E ^KTstl 

SC{l,...,Af} \ 



where 



(i + r^)(v/2io g (i + r' s 5) - >Ayr g ) 



(To = 4 

1 0, otherwise. 



(18) 



(19) 



Proof: Define random variables Y$ and Z5 as the 
following: 



1 



Y S = - log ^p) = 5 log(l+^s), for all SCM. 

(20) 

The facet defining constraints of C g (P,H) and C a (P) 

are of the form of 2~2ies Ri — an( l 2~2ieS ^ — 

E[Yg], respectively. Hence, by Definition [3] we have 
d H (C g (P, H),C a (P)) > 8 if and only if \Y S - E[Y S ]\ > 8, 
for all 5 C M. . After some manipulations, the following 
relations can be verified by employing Chebyshev's inequality: 



Pr 



{d H (C g (P,H),C a (P))>s} 



= Pr{max|y s -E[Fs]| > 



where a\ s denotes variance of Y s , and can be bounded from 
above as follows (c.f. Appendix II, [13]) 



^2 / u z s , 



l + Z s )U /21og(l + Z s 



(22) 



where 



J Z S 



(j2HiP^=r' s KT s . 



The desired result is concluded by substituting Z$ and u\ 
in ( 122b and combing the result with (fJTJi. ■ 

The system parameter <j\ in Lemma Q] is proportional to 
channel variations, and we expect it to vanish for small channel 
variations. The following lemma ensures that the distance 
between the optimal solutions of the utility maximization 
problem over two regions is small, provided that the regions 
are close to each other. 

Lemma 2: Let R\ and R* 2 be the optimal solution of 
maximizing the utility over C g (P,H±) and C g (P, H2), 



respectively. If there exist positive scalars A and B such that 

\u(Rx) - u(R 2 )\ < B\\R 1 -R 2 \\, 

for all R t e F(Cg(P, Hi)), i = 1,2. 

HR^-uiR.^ > AWRl-Rif, 

for all R t e C g (P, Hi), i = l,2, 



(23) 



and moreover if 

d H (C g (P,H 1 ),Cg(P,H 2 )) <8 
then, we have 



\R*-RZ\\<6? 



8* 



(*)' 



(24) 



Proof: Without loss of generality assume that u(R^) > 
u(R*). To simplify the notations for capacity regions, let C\ = 
Gg\P,H\) be a polymatroid, i.e., 



d = 1 R e R£ : Ri < f(S), for all S C M j>, (25) 

for some submodular function f(S), and let C2 be an expan- 
sion of Ci by 5 as defined in Definition |2] We first show that 
for every R E ^(Cz), there exists a vector R' E 3~{C\) such 
that ||U— R'\\ < 8, where T{-) denotes the dominant face of 
a capacity region as in Definition Q] 

Assume R is a vertex of C%- Then the polymatroid structure 
of C2 implies that R is the intersection of M constraints 
corresponding to a chain of subsets of M.. Hence, there is 
some k E M such that Rk = /({&}) + 8. Choose R' as 
follows 



(26) 



Ri — 8, i = k 
Ri , otherwise. 

R' is obviously in a (5-neighborhood of R. Moreover, the 
constraint corresponding to the set M. is active for R' , so we 
just need to show that R' is feasible in order to prove that it 
is on the dominant face. First, let us consider the sets S that 
contain fc. We have 

^i^=^^-5</(5). (27) 

ies ies 

Second, consider the case that k ^ S. 

R 'i = E + * 

i£SU{k} 

< f(SU{k})+8-R k 

< f(S) + f({k}) + 8-R l 
= f(S). 

where the first inequality come from (|2"7| >, and the second 
inequality is valid because of the submodularity of the function 
/(•)• 

The previous argument establishes that the claim is true 
for each vertex Rj of the dominant face. But every other 
point R on the dominant face can be represented as a convex 
combination of the vertices, i.e., 



E 



R = Z~l a 3 R ^ 



E a J 



l,a, > 0. 
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Using the convexity of the norm function, it is quite straight- 
forward to show that the desired R' is given by 

R ' = J2 a 3 R 'j> 



where Ftj is obtained for each Rj in the same manner as in 
d26]l. 

So we have verified that there exists a point, R, on the 
dominant face of C\ = C g (P, Hi) such that \\Rl - R\\ < S. 
By monotonicity of the utility function the optimal solution 
i?2 lies on the dominant face of C g (P, H2), hence, from the 
hypothesis and the fact that u(R%) 
conclude 



> u(Rl) > u(R), we 



u{R* 2 )~u(R) = \u(Rl)-u(R)\ < B\\Rl-R\\ < BS. (28) 

Now suppose that \\R\ — R\\ > (^S)i. By the hypothesis in 
( 1231 we can write 

u{Rl) -u(R) 



\u(R{)-u(R)\ > BS. 



(29) 



By subtracting 



from 



we obtain u(R%) < u(R*) 



which is a contradiction. Therefore, \\R\ — R\\ < (^S)i , and 
the desired result follows immediately by invoking the triangle 
inequality. ■ 
The following theorem combines the results of the above 
lemmas to obtain a bound on the performance difference of 
the greedy and optimal policy. 

Theorem 1: Let R* be the optimal solution to ([T3l . and 
TZ(-) the greedy rate allocation policy as defined in ([Pil l for 
a non-negative concave utility function Then for every 
e 6(0,1], 



u(JT) -u(E H [n(H)}) < eu(R*) + (1 - e)B(e) 



( OH 



B(e) 
A(e) 



( OH 



(30) 



where B(e) and A(e) are positive functions of e, such that for 
all H with d H (C g (P,H),C a (P))) < ^, they satisfy the 
following conditions. 



\u(Ra) - u(R g )\ 



\u(K(H)) - u(R) 



< 



> 



B(e)\\R a -R g \\, 

for all R a eT{C a {P)), 

for all R g € T{C g {P,H)), 

(31) 

for all Re C g {P,H). (32) 



Proof: Pick any e S (0, 1]. Define the event V as 

d H (C g (P,H),C a (P)))<^. 

By Lemma [T] the probability of this event is at least 1 
Conditioned on V, we have the following 



i(R*) -<u(K(H)) \ <B(e)\\TZ(H) - R*\\ 

- /<7H\k t (B(e) 



<B(e) 



A(e) 



( OH 



(33) 




0.2 



0.4 



0.6 



0.8 



Fig. 1 . Parameterized upperbound on performance difference between greedy 
and optimal policies. 



where the first inequality follows from monotonicity of the 
utility function and ( |3TT >. The second inequality is a direct 
result of applying Lemma [2] 

Using Jensen's inequality as in ( Tl3T > we can bound the left- 
hand side of d30l ) as follows 



u(R*)-u(E H [n(H)]) 

<u(R*)-E H [u(H(H))] 

< u(R*) - (1 - e)E H \u(iZ(H)) 



V 



-Pr{V c )E H u(U(H)) 



V 



< eu(R*) + (1 - e)E H \u{R*) - u{K(H))\ 



V . 

(34) 



In the above relations, the second inequality follows from 
Pr(V) > 1 — e, and the third inequality is obtained from 
non-negativity of the utility function after some manipulation. 
Replacing ( [331 in ( f34t gives the desired upperbound. ■ 

Theorem [T] provides a bound parameterized by e. For very 
small channel variations, gh tends to zero, and we can choose 
e proportional to an such that the bound in (f30b approaches 
zero. Figure Q] illustrates the behavior of the parameterized 
bound provided in d30b for different values of an- For each 
value of oh, the upperbound is minimized for a specific choice 
of e, which is illustrated as a dot in Figure [T] As demonstrated 
in the figure, for smaller channel variations tighter bound is 
achieved and the minimizer parameter decreases. 

The next theorem provides another bound demonstrating 
the impact of the structure of the utility function on the 
performance of the greedy policy. 

Theorem 2: Let R* be the optimal solution to ( fT3T > for the 
non-negative utility function u(R). Also let 1Z*(-) and TZ(-) 
be the optimal and greedy rate allocation policies, respectively. 
Then for every e <E (0, 1], 

u{R*)-u(E H [1Z(H)]) < eu{R*) + ^{l-e)r(e) 2 n, (35) 
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M 


H 




i=l 





where Q satisfies the following 

A max (-V 2 W (0) < n, for all £, ||£ - R*\\ < r(e), (36) 
and r(e) is given by 

r(e) = VI7^ + 
Ve 

i + E^i^ 

(37) 

Proof: Pick any e£ (0,1]. Define the event V similarly 
to the proof of Theorem Q] Because of monotonicity of the 
utility function, we know that R* lies on the dominant face 
of C a (P). Since the region C a (P) is the average of all regions 
C g (P, if), the optimal policy 1Z*(H) should give a point on 
the dominant face of C g (P, if), for almost all if. Therefore, 
conditioned on V, we can bound the set in which 1Z*(H) 
ranges, i.e., ||7?*(if ) — R*\\ < r(e), after some straightforward 
manipulations. Now let us write the Taylor expansion of the 
function u(-) at R*. We have 

u(R) = u(R*)+Vu(R*)'(R- R*) 

-±(R-R*)>(-V 2 u(t))(R-R*) 
> u(R*)+Vu(R*)'(R- R*) 

-±\\R- R*\\ 2 \ max (-V 2 u($)), 

for some £, ||£-ff*|| < ||ff-ff*||. 

(38) 

By replacing R by 72.* (if) and conditioning on V we have 
the following 

u{TZ*(H)) > u(R*) + 

Vu(R*)'(n*(H) - R*) - ^r(e) 2 n. 

Now we can bound the left-hand side of (f35T > by bounding 
the Jensen's difference u(R*) — Ejj [u(lZ* (if ))]. After some 
manipulation similar to ( f34t . we have 

u(R*) - u(E H [n(H)]) 

< u(R*)-¥. H [u(TZ*(H))] 



< u(R*) - (1 - e)E H u(K*(H)) 



V 



-Pr(V c )Eff u(Tl*{H)) 



V 



< u{R*) - (1 - e)(u(R*) - -r(e) 2 n 



eu{R*) + -{l-e)r{efn. 



Similarly to Theorem Q] Theorem [2] provides a bound 
parameterized by e which goes to zero for proper choice of e 
as Vl becomes smaller and the utility function tends to have 
a more linear structure. The behavior of this parameterized 
upperbound is also similar to the one illustrated in Figure [T] 

In summary, the performance difference between the greedy 
and the optimal policy is bounded from above by the minimum 
of the bounds provided by Theorem Q] and Theorem [2] 



V. Conclusion 

We addressed the problem of optimal resource allocation 
in a fading multiple access channel from an information 
theoretic point of view. We formulated the problem as a 
utility maximization problem for a more general class of utility 
functions. 

We considered two different scenarios. First, we assumed 
that the transmitters have power control feature and the 
channel statistics are known a priori. In this case, the optimal 
rate and power allocation policies are obtained by greedily 
maximizing a properly defined linear utility function. 

In the second scenario, power control and channel statistics 
are not available. In this case, the greedy policy is not optimal 
for nonlinear utility functions. However, we showed that its 
performance in terms of the utility is not arbitrarily worse 
compared to the optimal policy, by bounding their performance 
difference. The provided bound tends to zero as the channel 
variations become small or the utility function behaves more 
linearly. 

The greedy policy may itself be computationally expensive. 
A computationally efficient algorithm can be employed to allo- 
cate rates close to the ones allocated by the greedy policy. This 
algorithm just takes one iteration of the gradient projection 
method at each time slot. Under slow fading conditions, it 
can be shown that this method tracks the greedy policy very 
closely, and its performance is close to the optimal policy. 
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